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Abstract 

Using the theory of exponential Riordan arrays and orthogonal polynomials, we 
demonstrate that the "descending power" Eulerian polynomials, and their once shifted 
sequence, are moment sequences for simple families of orthogonal polynomials, which 
we characterize in terms of their three-term recurrence. We obtain the generating 
functions of the polynomial sequences in terms of continued fractions, and we also 
calculate their Hankel transforms. 
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1 Introduction 

The Eulerian polynomials [9, 14, 17, 21] 



PJx) 



k=0 



form the sequence Pn{x) which begins 

Po{x) = 1, Pi(x) = 1, P2{X) = 1 + X, P^ix) 

with the well-known triangle of Eulerian numbers [16] 
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as coefficient array. These coefficients Wn^k obey the recurrence [17] 

Wn,k = ik+ l)Wn-l,k + in- k)Wn-l,k-l 



1 



with appropriate boundary conditions. The closed form expression 

i=0 ^ ' 

holds. The polynomials Pjyx) were introduced by Euler [13] in the form 

fc=0 ^ ' 

They have exponential generating function 



71=0 



We have 

n 
fc=0 

and hence we can regard them as "descending power" Eulerian polynomials. 

In this note we show that the sequence of Eulerian polynomials Pni^x) is the moment 
sequence of a family of orthogonal polynomials. In addition, we show that the sequence 
of shifted Eulerian polynomials P„+i(x) is similarly the moment sequence of a family of 
orthogonal polynomials. For this, we will require three results from the theory of exponential 
Riordan arrays (see Appendix for an introduction to exponential Riordan arrays). These are 
[5, 6] 

1. The inverse of an exponential Riordan array [(?,/] is the coefficient array of a family 
of orthogonal polynomials if and only if the production matrix of [(7, /] is tri-diagonal; 

2. If the production matrix of [(7, /] is tri-diagonal, then the elements of the first column 
of [(?, /] are the moments of the corresponding family of orthogonal polynomials; 

3. The bivariate generating function of the production matrix of [(7, /] is given by 

e'y^Zix) + A[x)y) 

where 

A{x) = /'(/»), 

and 

where f{x) is the compositional inverse (series reversion) of f{x). 
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A quick introduction to exponential Riordan arrays can be found in the Appendix to this 
note. For general information on orthogonal polynomials and moments, see [8, 15, 28]. 
Continued fractions will be referred to in the sequel; [29] is a general reference, while [18, 19] 
discuss the connection of continued fractions to orthogonal polynomials, moments and Hankel 
transforms [20, 25]. We recall that for a given sequence a„ its Hankel transform is the 
sequence of determinants hn = |aj+j |o<jj<n- Many interesting examples of number triangles, 
including exponential Riordan arrays, can be found in Neil Sloane's On-Line Encyclopedia of 
Integer Sequences [26, 27]. Sequences are frequently referred to by their OEIS number. For 
instance, the binomial matrix (Pascal's triangle) B with (n, /c)-th element is A007318 . 



2 The Eulerian polynomials Pn{x) 



We consider the sequence with e.g.f. 



(a - /3)e("-'3)* 



This is the sequence that begins 



1, a, a{a + a{a^ + 4a/3 + a{a^ + iWp + lla/3^ + . . . . 



Setting a = 1 and (3 = x gives us the Eulerian polynomials P„(x). We have the 



Proposition 1. The production matrix of the exponential Riordan array 



{a - /3)e("-^)* e("-^)* - 1 



is tri-diagonal. 



Proof. Writing the above exponential Riordan array as [g, /], we have 



fit) 



g(Q-/3)t _ ]_ 



and hence 



nt) 



and 




Then 



A{t) = f'{f{t)) = (at + l)(/3t + 1) = 1 + (« + /3)t + ape. 



We have 



git) 



{a - /3)e("-^)* 
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and hence 



and so 



g'it) 



gim) 

Thus the production matrix, which has bivariate g.f. given by 

e^y{a + apt +(! + (« + + al3t^)y), 

is tri-diagonaL 

We note that the production matrix takes the form 



□ 



( a I 

a/3 2a + /3 1 

4a/3 3a + 2/3 1 

9a/3 4a + 3/3 












1 

16a/3 5a + 4/3 1 

25a/3 6a + 5/3 



v ■■ 
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For completeness, we note that while in the special case a = /3 the Riordan array is not 
obviously well-defined, the production matrix is, and it leads in this special case to the 
exponential Riordan array 

1 t 



1 — at 1 — at 

which has general element (^) Ija"""^. In the case a = /3 = 1, we get the exponential Riordan 
array 



-t l-t 

whose inverse is the coefficient array of the Laguerre polynomials [3]. 

Returning now to the Eulerian polynomials, we set a = 1 and /3 = x, to get 

Theorem 2. The Eulerian polynomials Pn{x) are the moments of the family of orthogonal 
polynomials Qn{t) defined by Qo(t) = 1, Qi{t) = t — 1, and 

Quit) = {t- {{n - l)x + n))g„_i(t) -in- lfxQn-2{t). 

Proof. The initial polynomial terms of the sequence Qn{t) can be read from the elements of 



xe 



[l-x)t ' ]_ _ ^g{l-x)i 



l + t I 



In 



X 



l + t 
l + xt 
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which begins 



/ 1 

-11 

2 -x-3 1 

-6 2a;2 + 5x + ll -3{x + 2) 1 

V ; ; ; ; 



Hence in particular Qo{t) = 1 and Qi{t) = t — 1. The three-term recurrence is derived from 
the production matrix, which in this case is 
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4 + 3x 


1 
















16x 


5 + 4x 
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25x 


6 + 5x . 



















□ 

Corollary 3. The sequence of Eulerian polynomials Pn{x) has ordinary generating function 
given by the continued fraction 

1 



1-t 



xr 



1 - (2 + x)t 



AxV 



1 - (3 + 2x)t - 



9xt' 



Corollary 4. The Hankel transform of the sequence of Eulerian polynomials Pn{x) is given 
by 



h^ = xr.)i[k\\ 



k=l 



3 The shifted Eulerian polynomials Pn+i{x) 

For the shifted Eulerian polynomials P„_|_i(x), we consider the exponential Riordan array 

[9'it)J{t)], 

where 



{a - /3)2e(°+^)* 



where we retain the use of g{t) = ^'^Z^J^ia-m't from the previous section 



(a-/3)e("-/3)* 

When a = 1 and /3 = x, g'{t) generates the shifted sequence P„+i(x). We then have 
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Proposition 5. The production matrix of the exponential Riordan array 
is tri- diagonal. 

Proof. As in the previous proposition, we obtain 

Mt) = f'ifit)) = + + 1) = 1 + (a + /3)t + al5t\ 

where 



m = 



In 



at + 1 



Then 



Z{t) 



a-P \/3t+l 

g'Vit)) 



a + /3) + 2al3t. 



9'im) 

The bivariate generating function of the production matrix is then 

e*^((a + + 2a/3t + (! + (« + (3)t + a(3t^)y), 
and hence the production matrix is tri-diagonal. 
The production matrix in this case begins 



/ a + 13 1 

2a/3 2(a + /3) 1 

6a/3 3(a + /3) 1 

12a/3 4(a + /3) 

V : ! : : 

In the case a = /9, we obtain the exponential Riordan array 

1 t 



\ 



with (n, A;)-th element (^+J) f 



n+l\ n!Q,n fc_ p^j, q, = ^ = ^ ^]^ig gives US 



□ 



which is A105278 . 

Specializing to the values a = 1, /3 = x, we get the 

Theorem 6. The shifted Eulerian polynomials Pn+i{x) are the moments of the family of 
orthogonal polynomials Rn(t) given by Ro{t) = 1, Ri{t) = t — x — 1, and for n > 1, 

Rn{t) = {t- n{l + x))Rn-iit) - n{n - l)xi?„_2(t). 
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Proof. The initial terms of the polynomial sequence Rn{t) can be read from the elements of 
the inverse matrix 



(a - /3)2e(-+/5)* e(i-")* - 1 



-, -1 



In 



[l + t){l + tx)' 1- X \l + xt 



1 + t 



which begins 



1 

-X - 1 1 

2x2 + 2a; + 2 -3(x + 1) 1 

-6{x^ + x"^ + X + 1) llx^ + Ux + 11 -6(x + l) 1 
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The three-term recurrence is derived from the production matrix, which in this case is 
/l + x 1 ...\ 



+ X 


1 














2x 


2(1 + x) 
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3(1 + x) 
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6(1 + x 



□ 

Corollary 7. The sequence of shifted Eulerian polynomials P„+i(x) has ordinary generating 
function given by the continued fraction 

1 



l-(l+x)t 



2xt' 



1 -2(1 



6xt' 



l-3(l + x)t 



12xt' 
1 - ■ ■ 



Corollary 8. The Hankel transform of the shifted Eulerian polynomials Pn+i{x) is given by 



-n- I k + 1 

2 



n—k 



/.„ = (2x)mn 

fe=l 

4 The Eulerian number triangles 

As with the Narayana numbers and their associated number triangles [7], we can distinguish 
between three distinct but related triangles of Eulerian numbers. Thus we have the triangle 
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A173018 [16, 17] 

/ 1 ... \ 

1 ... 

1 1 ... 

1 4 1 ... 

1 11 11 1 ... 

1 26 66 26 1 ... 

\\ \ \ \ \ \ ■■. j 
of Eulerian numbers Wn^k that obey the recurrence 

Wn,k = {k + l)Wn-i,k + {n- k)Wn-l,k-l 
with appropriate boundary conditions, for which the closed form expression 

n-k , _|_ 1 \ 

i=0 ^ ' 

holds. We have the reversal of this triangle, which is the triangle A123125 of the coefficients 
Ar,^k [1] where 

An.k = Y^{-\y["'^.\h-^)\ 

which begins 
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and finally we have the Pascal-like triangle of coefficients 
which begins 
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This is A008292 . 
We have 

An,k = {n-k + l)An^ik - 1 + {k + l)i„_i/c, 

with appropriate boundary conditions. As with the Narayana numbers, each of these tri- 
angles has signficant combinatorial applications and it is often important to distinguish one 
from the other. 

Example 9. The sequence a„ = ^^=o ^n,k^^ is the sequence A000670 of preferential ar- 
rangements, or rankings of competitors in a race, with ties [23]. The sequence 



k=0 k=0 



or A000629 is the sequence of rankings of competitors in a race, with ties and dropouts [22]. 
Note that from our results above, the sequence has generating function given by 



2x' 

1 — X 



1-Ax 



Ix^ 



18x2 

1 - 7x 



1 - 

The g.f. of the sequence a„+i is given by 

1 



Ax' 

1 - 3x 



I2x 

1 — 6x 



24x2 
1 - 9x 



In this case it happens that 6„ is the binomial transform of a„, and hence [4] its g.f. has 
continued fraction expression 

1 



2x' 

l-2x 



1 — 5x 



Ix^ 



18x2 

1 - 8x 



1 

5 A related ODE 

The form of f{t) above is related to a simple ODE. This arises as follows. In order to have 
a tri-diagonal production matrix, we need to have an expression of the form 
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Now substituting z 



fit) we obtain 



f{f{fm = i+fifit) + uf{tf 



or 



f(t) = l + fif{t) + uf{tf 



or 



— = l + fiy + uy , 
where y = f{t). In the Eulerian case above, we have 

^ = (l + «y)(l + /?!/), 



with initial condition |/(0) = 0. The form oi y = f(t) follows from this variant of the logistic 
equation. 

6 Appendix: exponential Riordan array 

The exponential Riordan group [2, 10, 12], is a set of infinite lower-triangular integer matrices, 
where each matrix is defined by a pair of generating functions g{x) = go + gix + g2x'^ + • • • 
and f{x) = fix + /2X^ + • ■ ■ where go 7^ and /i 7^ 0. We usually assume that 



The associated matrix is the matrix whose i-th column has exponential generating function 
g{x)f{xy/i\ (the first column being indexed by 0). The matrix corresponding to the pair 
/, g is denoted by [g, /]. The group law is given by 



The identity for this law is / = [1, x] and the inverse of [g, f] is [g, f]~^ = [1 / {g o f) ^ f] where 
/ is the compositional inverse of /. 

If M is the matrix [g, /], and u = (u„)„>o is an integer sequence with exponential generat- 
ing function U (x), then the sequence Mu has exponential generating function g[x)U{f{x)). 
Thus the row sums of the array [g^ /] have exponential generating function given by g{x)e^^^^ 
since the sequence 1,1,1,... has exponential generating function e^. 

As an element of the group of exponential Riordan arrays, the binomial matrix B with 
(n, /c)-th element (^) is given by B = [e^,x]. By the above, the exponential generating 
function of its row sums is given by e^e^ = e^^, as expected (e^^ is the e.g.f. of 2"). 

To each exponential Riordan array L = [g, f] is associated [11, 12] a matrix P called its 
production matrix, which has bivariate g.f. given by 



90 = fi = 1- 



[gJ]-[h,l] = [gihof)Jof]. 



e^^(Z(a;) + A{x)y) 



where 



= /'(/»), Z{x) 



9im) ■ 
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We have 

P = L-^L 

where L [24, 29] is the matrix L with its top row removed. 
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